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Neutrino helicity asymmetries in leptogenesis 
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It is pointed out that the heavy singlet neutrinos characteristic of leptogenesis develop asymmetries 
in the abundances of the two helicity states as a result of the same mechanism that generates 
asymmetries in the standard lepton sector. Neutrinos and standard leptons interchange asymmetries 
in collisions with each other. It is shown that an appropriate quantum number, B — L' , combining 
baryon, lepton and neutrino asymmetries, is not violated as fast as the standard B — L. This 
suppresses the washout effects relevant for the derivation of the final baryon asymmetry. One 
presents detailed calculations for the period of neutrino thermal production in the framework of the 
singlet seesaw mechanism. 
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I. INTRODUCTION 



Leptogenesis is an attractive way of generating the 
baryon number of the Universe 0. The main idea put 
forward by Fukugita and Yanagida 0, 0, is that a 
lepton number asymmetry can be produced in the de- 
cays of heavy singlet neutrinos into leptons and Higgs 
bosons and such an asymmetry is partially transferred 
to the baryon sector through electroweak sphaleron pro- 
cesses 5] that violate B and L but not B — L. The 
mechanism requires nonconservation of lepton number 
and CP provided by neutrino Majorana masses and com- 
plex Yukawa couplings. Both masses and couplings form 
the well known singlet seesaw model 

a I la E3 of 

light neutrino masses and thus establish a close relation- 
ship between baryogenesis and low e nergy p henomenol- 
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l25l |26| . This connection contributed for the present wide 
interest in leptogenesis. 

The calculation of the final baryon asymmetry has 
been done in the literature with increasing levels of ac- 
curacy [II [II il il il El IS |H H3 but the main 
elements have remained the following. Singlet neutrino 
reactions are not symmetric under CP due to nontrivial 
complex Yukawa couplings in the neutrino mass eigen- 
state basis. Departure of neutrino densities from thermal 
equilibrium values are a necessary condition 31] to ob- 
tain net lepton asymmetry sources. This occurs if neu- 
trinos are not produced in the inflaton decay but only 
gradually from active lepton and Higgs boson collisions. 
It occurs also to some extent when any of the neutrino 
species undergoes the transition to the respective nonrel- 
ativistic temperature epoch. Weak sphaleron processes 
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transform a fraction of the generated B — L asymme- 
try into baryon number. The final B — L and baryon 
asymmetries depend on lepton number violating reac- 
tions whose net effect is to dissipate B — L. They include 
the reactions 4>4> — ► Z,Zj, <plj — > tpU, hlj — > <p4>, that vio- 
late lepton number by two units, but also processes that 
violate standard lepton number by one unit [3j such as 
top quark and electroweak gauge boson scatterings like 
tqt -> N a U, 4>W -> N a k and (f> B -> N a k. Neutrino 
number densities and the set of standard lepton, quark 
and Higgs boson number asymmetries obey a system of 
coupled Boltzmann equations that is necessary to inte- 
grate to derive the final (present) B — L asymmetry. 

And here comes the main point of this paper. This 
system is incomplete because, contrary to what has ever 
been assumed, the two helicity states of singlet Majo- 
rana neutrinos do not have exactly the same abundances. 
We make the case that leptogenesis mechanisms natu- 
rally generate asymmetries in the two neutrino helicity 
state abundances, in the same way as the standard lep- 
ton asymmetries, and that the neutrino helicity asym- 
metries play a role in the system of Boltzmann equations 
that govern the evolution and transport of particle asym- 
metries, lepton number in particular, and therefore con- 
tribute to the determination of the final B — L asymmetry 
as a function of the fundamental parameters of the the- 
ory. We present explicit calculations in the framework of 
the singlet seesaw mechanism. 

In next section we review some of the properties of 
this model that are relevant for lepto gen esis in relation 
with the light neutrino mas s sp ect ra 13^. l33l l34| as in- 
dicated by solar |3lllMl3^l3al39ll40n4l|. atmospheric 
|H OH El an d terrestrial neutrino experiments 0, ^(| . 
In section llHl we introduce the concept of neutrino helic- 
ity asymmetries and discuss their relevance for the trans- 
port of lepton number. An appropriate quantum number 
B — V is proposed to replace the usual B — L difference. 
In section llVl we calculate explicitly the lepton and neu- 
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trino asymmetries generated during the phase of singlet 
neutrino thermal production. This was presented in a 
brief fashion in a meeting [i^ . Here we give a complete 
account of the work and improve the integrations over 
phase space by including Pauli blocking and lepton ther- 
mal mass effects in the numerical calculations. The Higgs 
boson thermal mass had already been taken into account. 
In section we study the washout processes and evalu- 
ate the damping rate of B — L'. We compare with the 
traditional treatment without neutrino helicity asymme- 
tries and take the appropriate lessons. The main results 
are summarized in the final section. 



II. SEESAW MODEL 

The singlet seesaw mechanism [(J adds to the standard 
model singlet (left-handed) neutrinos, N a , with heavy 
Majorana masses and Yukawa couplings with the stan- 
dard lepton and Higgs doublets, U and <f>, of the form 



hi a kN a (f) + -M a N a N a + H.C. 



(1) 



Spontaneous breaking of SU(2) x U(l) yields the light 
neutrino mass matrix (v = (</> )) 



-(hM- 1 h i ) ij v 



(2) 



The proper decay rate of N a into leptons and Higgs is 
given by 



1 N a 



8tt 



-Mn 



(3) 



if one ignores thermal effects. Delayed decay occurs when 
the ratio to the Hubble expansion rate H at the temper- 
ature T = M a , 



K„ 



1 A 



N„ 



H, 



(T=M a ) 



(4) 



is small. In the radiation era, H = 1.66 g\ T 2 /Mp, 
where denotes the number of relativistic degrees of 
freedom, 107.5 in the standard model. It is enough to 
compare the sum 



K = Y^ K a = (10 3 cV" 1 ) Tr^M" 1 /^ 2 



(5) 



with the light neutrino mass scale Tr[m] to conclude that 
the delayed decay condition, K a < 1, is in general in 
conflict with the atmospheric neutrino mass gap |42| . 
Am 2 w 2.5 ■ 10~ 3 eV 2 , which implies K > 50. Strictly 
speaking, the delayed decay scenario only requires that 
the lightest of the heavy neutrinos satisfies K a < 1. But 
that is not the most natural picture, in particular if light 
neutrinos are quasidegenerate. In section [Tvl we will as- 
sume that all parameters K a are large, of the order of 
50 or more, which has the effect that singlet neutrinos 
enter in thermal equilibrium at relativistic temperatures 

T a » M a . 



III. NEUTRINO HELICITY ASYMMETRIES 

Sterile neutrino Yukawa couplings UN a (j) conserve total 
lepton number, Lt, assigned as L = — 1 for left-handed 
neutrino fields, N a , and L — +1 for the right-handed 
neutrino conjugate fields N a . Neutrino Majorana masses 
break lepton number but it is clear that in the ultrarela- 
tivistic limit the masses are negligible and neutrinos and 
antincutrinos may have unequal abundances and sym- 
metric chemical potentials like any other particles. In 
that case neutrinos carry lepton number that can be ex- 
changed between them and standard leptons in collisions 
mediated by Yukawa interactions. The Majorana masses 
do not change this completely. One faces a similar prob- 
lem at defining lepton number of Majorana mass solar or 
atmospheric neutrinos, or lepton asymmetries of neutri- 
nos at the time of Big Bang Nucleosynthesis |48|. The 
main difference is that they deal with standard active 
neutrinos and not necessarily sterile neutrinos. Free mas- 
sive neutrino states are solutions of the Dirac equation 
that can be discriminated as spin eigenstates and in par- 
ticular as helicity eigenstates. The ±1/2 helicity eigen- 
states spinors satisfy the relation 



1 / E-Mj" 
- 1 ±75 

2 V V 



(6) 



where E is the neutrino energy, p the momentum and M 
the mass. When neutrinos are ultrarelativistic, helicity 
and chiral states are almost identical and their lepton 
number is maximal (±1). For an arbitrary neutrino state 
the lepton number current density is evaluated as the 
expectation value 



(7) 



where N is either a left-handed chiral field N a or 
its conjugate N a , and x = Y a + N a is a Majorana 
field (for the quantization of Majorana mass fields see 
refs. HE HIJ Ull)- As a result a positive (negative) helicity 
eigenstate carries a well defined average lepton number 
equal to the neutrino speed v — p/E (— v): 



L = — ; — = ±v 
u'u 



(8) 



The lepton number vanishes in the nonrelativistic limit 
due to the Majorana nature of the neutrino mass. 

One may observe that in contrast to lepton number, 
helicity is not invariant under Lorentz transformations. 
However, in an isotropic Universe the comoving thermal 
bath frame is a privileged frame where isotropy enforces 
the spin density matrix to be diagonal in the helicity 
basis. That means that each of the neutrino flavors A^ a 
can be divided in two populations of opposite helicities 
and well defined distribution functions /„ . The total 
lepton number carried by each neutrino species is equal 
to 

V I -/-)«„, (9) 
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where v a is the neutrino speed and V the spatial volume. 
It will prove convenient to work instead with the helic- 
ity asymmetries i.e., the differences between positive and 
negative helicity neutrino abundances, 



A„ 



V 



d 3 Pa 

(2tt) 3 



(10) 



Sterile neutrinos have been deprived of lepton number 
in the leptogenesis literature because they have Majo- 
rana masses (one exception is the oscillation mechanism 
of ref. |52(). But contrary to the lepton number assign- 
ment, that is to an extent arbitrary, it is a unambiguous 
fact that neutrinos develop helicity asymmetries as a re- 
sult of collisions with standard leptons and Higgs and 
also directly from the same leptogenesis processes that 
generate lepton asymmetries. That is shown explicitly in 
section HvTI 

It is important to realize that neutrino helicity asym- 
metries affect the way lepton number is transported in 
collision processes and consequently the so-called lepton 
number washout effect. Take for example scatterings like 
N a li <-> iq tl N a k <-> tqt, and crossed channels, that vio- 
late standard lepton number by one unit (t is the right- 
handed top quark and qt the quark iso-doublet linear 
combination that has a Yukawa coupling with it). The 
correct evaluation of the lepton number violation rate 
must take into account that the reaction rates depend on 
the neutrino helicity and the two helicity states have dif- 
ferent abundances. Reactions with positive helicity neu- 
trinos in initial or final states like iV+Z, are suppressed 
with respect to the opposite helicity states N~U. The for- 
mer are possible only because neutrinos have Majorana 
masses and their transition amplitudes are suppressed by 
the Lorentz contraction factors M a /E. If neutrinos are 
in thermal equilibrium their helicity asymmetries can be 
parametrized with degeneracy parameters rj a = fJ, a /T as 
any other fermions, r\ a for positive helicity and — r\ a for 
negative helicity states, and the rate of standard lepton 
number violation contains terms like 

lN a U->tq t ~ Itqt^Nah = J+(V<f> + Vh + Va) 

+ l-iVt + Vli ~Va) , (11) 



where 7+, 7_ stand for the average reaction rates of 
N+k <-» tq t and N~U <-* iq t respectively. We have re- 
placed rjt and rj qt with the Higgs degeneracy parameter 
using the constraint rj t — rj qt = enforced by the rapid 
top quark Yukawa interactions. Identical expressions ap- 
ply to other single neutrino absorption and emission re- 
actions, decays and inverse decays N a <-> li<p (at tem- 
peratures T < M a where the Higgs and lepton thermal 
masses are small enough), or <f> «-» N a li (at temperatures 
T > 2M a ), and respective radiative processes with one 
electroweak gauge boson in the initial or final state. No- 
tice that because 7+ is smaller than 7_ the result depends 
on the neutrino chemical potentials. The two rates co- 
incide with each other only when neutrinos are at rest. 
On the other hand one should not expect that the neu- 
trino degeneracy parameters rj a are damped faster by the 



lepton number violating reactions than the combinations 
V4> + Vh 1 a l s o damped by AL — 2 collisions like cf>li <j)lj . 

Contrary to the partial lepton numbers of standard 
leptons, Li, and neutrinos, L a , the total lepton number 



L'i 



(12) 



and the quantum number B — Lt are conserved by 
neutrino- lepton Yukawa couplings. In addition, B — Lt is 
conserved by sphalerons and is only violated by neutrino 
Majorana masses. For practical purposes namely, ap- 
plication of the constraints imposed by CPT invariance 
and unitarity on the collision rates, it is more convenient 
to work with the helicity asymmetries A a rather than 
the lepton numbers L a because the states with definite 
distribution functions are helicity eigenstates, not (chi- 
ral) lepton number eigenstates. On the other hand, the 
quantum number 



(13) 



has essentially the same interesting properties as the to- 
tal lepton number Lt- The combination B — L', also 
conserved by weak sphalerons, reduces to the standard 
model B — L when the heavy neutrinos vanish from the 
Universe. It is conserved by neutrino Majorana masses 
but its violation by Yukawa couplings is suppressed by 
the neutrino masses i.e., mass over energy ratios. 

Neutrino helicity asymmetries play a role in the trans- 
port of standard lepton number. We have shown for in- 
stance that processes that violate flavor quantum num- 
bers Li and A a by one unit contribute to the violation 
rate of Li — B/3 as follows (if all particles are in thermal 
equilibrium) : 



d{Lj - g/3) 
dt 



-lfa\v4> 
-l?a\v<P 



Vh 
Vh 



Va) 
Va) 



(14) 



where 7^ is the total rate of L' conserving reactions 

and 7 la the total rate of AL' = ±2 reactions. In turn, 
the standard lepton asymmetries are transferred to the 
neutrino sector in collisions with them. The same pro- 
cesses as above contribute to the violation rates of helicity 
asymmetries A a as 



dA a 
dt 



(V<t>+Vh~Va)- ifa iV<t> + Vl i + Va ) + " ' 



(15) 



This proves that the neutrino helicity asymmetries can- 
not be assumed identically zero because the L' conserving 
reactions are faster than the L' violating ones. It is also 
clear that only the later contribute to the violation rate 
of B — L' . We will return to this point in section IVI 

This discussion shows that in order to correctly eval- 
uate the lepton number washout effects neutrino helic- 
ity asymmetries are an essential ingredient. They have 
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been so far completely ignored in the leptogenesis liter- 
ature. The calculation of the lepton number generated 
during the decay phase of the lightest neutrino (s) and in- 
tegration of the Boltzmann equations including neutrino 
hclicity asymmetries is complicated by the fact that the 
neutrinos are neither purely nonrelativistic nor ultrarel- 
ativistic at temperatures close to their masses. In this 
paper we limit ourselves to the temperature range where 
neutrinos are ultrarelativistic which permits first order 
calculations on their mass over temperature ratios. In 
next section we calculate the B — L' asymmetry gener- 
ated during neutrino thermal production. 

IV. GENERATION OF LEPTON 
ASYMMETRIES 

A. Neutrino thermal production 

Let us now examine the leptogenesis processes in de- 
tail during neutrino thermal production. We assume that 
the Universe is initially empty of singlet neutrinos, not 
produced in the inflaton decay but only thermally from 
standard leptons and Higgs. The dominant thermalizing 
reactions are identical to the ones of the charged lep- 
ton isosinglets |53|. namely, top quark scattering pro- 
cesses like qtin — ► i,iV a , Higgs boson decay (j> — > kN a 
and related scattering processes with one additional elec- 
troweak gauge boson in the initial or final state. Higgs 
boson decays into leptons and neutrinos, first considered 
in ref. [42] in the context of leptogenesis, are allowed 
as much as the decays into a lepton iso-doublet and a 
charged lepton isosinglet, <f> — > ^ej, because the Higgs 
has a significant thermal mass, rn& = x^T ~ 0.6 T, larger 
than the lepton thermal masses |5J,|55|,|56| as pointed out 
in ref. 0. 

Let n a denote the average neutrino number densities 
per helicity degree of freedom and Y a = n a V the abun- 
dances in a fixed comoving volume whose spatial vol- 
ume V expands as V = 3HV. For definiteness we as- 
sume that the reheating temperature is much higher than 
the neutrino masses so that neutrinos thermalize while 
they are ultrarelativistic. Then, the equilibrium densi- 
ties and abundances are equal to n oq = 0.90 T 3 /n 2 and 
Y C q — n cq V, respectively. Assuming for simplicity that 
the distribution functions scale with the equilibrium dis- 
tribution functions /„ q as 

77 

fa = — /r , (16) 

and approximating the Pauli blocking factors 1 — f a with 
1 — fa q in neutrino emission reactions, the neutrino abun- 
dances evolve as 

Y a = T a (Y c(i -Y a ) . (17) 

Here we neglect reactions with two or more neutrino 
states that are of higher order in the Yukawa couplings. 



T a represent the neutrino collision frequencies. At tem- 
peratures much higher than the neutrino masses T a scale 
with the temperature as follows: 

r a = -?-(h)h) aa T. (is) 

They are comparable with the <f> — > hN a proper de- 
cay rate, equal to (/it/i^m^ if one ignores ther- 
mal effects other than the Higgs boson mass. The coef- 
ficient [3 gets contributions from the Higgs decay, top 
quark and W, B gauge boson scatterings in analogy 
with the charged leptons case [2^, [5j|. The relative 
weights of these reactions depend on the temperature be- 
cause the coupl ings constants |57| and thermal mass fac- 
tors 0, IHBl run with the energy scale J25| . Using 
— 0.6 T for the Higgs mass and 10 9 GeV temperature 
scale coupling constants namely, a s ~ 1/26, a w « 1/38, 
a' 1/81 for the strong, SU(2) and U(l)y gauge inter- 
actions respectively, and A t w 0.60 for the top Yukawa 
coupling, one obtains f3 w 1/7, where 46% comes from 
the Higgs boson decay, 41% from W and B electroweak 
gauge boson scatterings, and 13% from right-handed top 
quark scatterings. 

In the radiation era the Hubble expansion rate scales as 
H = l/2t oc T 2 and the assumed initially zero neutrino 
densities n a converge exponentially to the equilibrium 
densities: 

n a = n oq (l - e~ T «/ T ) . (19) 
The relaxation temperatures T a are given as 

T a = ^ = [3K a M a , (20) 
ti 

and the second identity establishes a relation with the 
parameters K a of Eq. 10 that control the speed of neu- 
trino decays when they become nonrelativistic and vanish 
from the Universe. The low energy neutrino data indi- 
cates that the parameters K a > 50 are large, see Eq. |JSJ . 
As a result singlet neutrinos reach thermal equilibrium 
at temperatures T a > 10AT a when they are still ultrarel- 
ativistic. 

We have defined neutrino densities in the basis of Ma- 
jorana mass eigenstates as usual. However, this is not 
valid for all temperature scales because neutrinos also 
get thermal masses from the interactions with the lep- 
ton - Higgs thermal bath. The chiral conserving thermal 
mass terms are [3 E3 

m\ h = \{tih) ab T 2 . (21) 

At high enough temperatures the thermal masses domi- 
nate over the vacuum masses M a and the neutrino Hamil- 
tonian eigenstates are eigenstates of the matrix (h^h) a b- 
The opposite happens when M a are much larger than 
the thermal masses, and at temperatures where the two 
types of masses are comparable with each other the neu- 
trinos undergo strong flavor oscillation processes. One 
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can estimate the thermal over vacuum mass ratios using 
Eqs. 0, (Hi: 



K\M a T 
10 19 GeV T 2 



(22) 



These have to be small at the relaxation temperatures T a 
in order that the neutrino densities evolve as described 
in this section. For decay constants K a of the order of 
50 (100) this happens for neutrino masses M a below 10 12 
(10 11 ) GeV. Then, the thermal masses can be neglected 
provided that the vacuum masses are not degenerate. 

There is another point. The neutrino states produced 
in collisions are linear combinations of the mass eigen- 
states N a of the form hi a N a . In a free path the mass 
eigenstate wave functions oscillate with unequal frequen- 
cies, -\/p 2 + M 2 , for the same linear momentum p. They 
differ by (M 2 — M 2 )/2_E for relativistic particles (energy 
E w p) and in a mean free path 1 /Y a give rise to average 
phase differences equal to 



M, 



Ml 



2T a 



Ml 



Ml 



(23) 



where we made use of Eqs. Il2()f and with (3 as 1/7. 
This shows that as long as thermal masses are much 
smaller than the vacuum mass gaps the neutrino mass 
eigenstates develop large phase decoherence between col- 
lisions so that one can consider that the neutrino states 
with definite number densities are the vacuum mass 
eigenstates N a . This is important because if, for example, 
the neutrino densities n a were all equal to each other the 
leptogenesis sources would vanish in the ultrarelativistic 
regime at lowest order, as the results of next sections 
show. 



B. Standard lepton asymmetries 

Leptogenesis is dominated by the following CP asym- 
metric reactions j^: Higgs decays into leptons and sin- 
glet neutrinos, inverse decays and scatterings of leptons 







FIG. 1: Diagrams contributing to CP-asymmetries in decays, 
inverse decays and scatterings. 



off neutrinos. The CP asymmetries result from the di- 
agrams of Fig. more specifically, from the interfer- 
ence between the tree level amplitude and the absorptive 
part of the one-loop amplitude of the Higgs boson decay 
(inverse decay) and from the interference between the t 
channel amplitude and the absorptive part of the s chan- 
nel amplitude of neutrino-lcpton scattering. The particle 
asymmetry sources vanish in thermal equilibrium but, 
as long as neutrinos stay rarefied scatterings and inverse 
decays do not match Higgs decays and particle asymme- 
tries develop in the various lepton flavors, singlet neutri- 
nos, and Higgs boson as well as enforced by hypercharge 
conservation. 

The source terms (labeled with S) responsible for the 
generation of the standard family lepton numbers Li are 
given in leading order by 

(Li) s = A 7 (0 -> N b k) s ~ A 7 (N b k -> $s 

b 

+ J2 A l( N Jj -» N b k)s ~ Aj(N b U -» NJ^s , (24) 
abj 



where 



A 1 {X^Y)= 1 {X^Y)- 1 {X^Y) (25) 



denotes the difference between the rate of an arbitrary 
reaction X — > Y integrated over a fixed comoving vol- 
ume and the rate of the conjugate reaction X — > Y, 
where X and Y are CPT conjugate states of X and 
Y respectively (CPT transforms particles into antipar- 
ticles, reverses helicity but not momentum). The dif- 
ference Aj(X — » Y) can be separated in two types of 
contributions, the transport terms proportional to par- 
ticle antiparticle abundance asymmetries and the source 
terms responsible for primordial asymmetry generation. 
By definition the source terms exist in absence of particle 
aymmetries and in the Aj(X — > Y)s expressions we as- 
sume that CPT conjugate states have exactly the same 
distribution functions. 

CPT invariance and unitarity ensure that the total de- 
cay rates of any quantum state and its CPT conjugate 
state are equal to each other. This translates in con- 
straints on the rate asymmetries: 



(26) 



In leading order, one gets: 

A 7 (N b k -> 0) s + &l(N b k -» NJ^s = , (27) 



J2 A 7 (0 -» N b h) s = 



(28) 



Then, the Li source terms of Eq. (|24|l can be written as, 
(Li)s =J2A^^N b l i ) s + J2 A 7(N a l j ^N b l i )s.(29) 

b abj 



6 



On the other hand, CPT invariance implies that any 
reaction X — > Y has the same transition probability as 
Y — > X where the bars indicate CPT conjugate states. 
As a result the rate asymmetries A 7 (A — > y) obey 
the following constraints in thermal equilibrium assuming 
identical particle and antiparticlc distribution functions 
(zero chemical potentials): 

a 7 (a y) cq + A 7 (y A) cq = o . (30) 

Here one employs the identity jES E3 satisfied by the 
thermal distribution functions of particles involved in two 
inverse reactions X — > Y and Y — > X, schematically, 
/x(l ± /y) = ± /jc), where f x (fy) stands for the 
product of initial state particle distribution functions and 
1 ± fy (1 ± fx) for the product of final state stimulated 
emission and/or Pauli blocking factors. Applying the 
above constraint on Eqs. I|27|) and l|28[l one obtains 

A 7 (0 JV 6 Zi)eq + A 7(^ ~> Wi)eq = 0, (31) 

^(N b k 0) oq = ^ A 7 (JV„f J - - = 0. (32) 

bi abij 

This ensures that Eq. (|29H satisfies the well known prop- 
erty that the as ymm etry source terms vanish in thermal 
equilibrium 59, 60j. 

The diagrams of Fig. JIJ lead to the following results: 

A 7 (0 -» N b k)s = 

£jy«6 / d * UmFaF b (p -_^ 2 P l ml , (33) 

A 7 (AT Q Z j _> jv&) s = 

Jijab J d<P fafjF^FiF b (pa ^p_ m , , (34) 

where f a are particle distribution functions, F a are 
fermion Pauli blocking factors 1 — f a or the stimulated 
emission factor i 7 ^ = 1 + and 

is the phase space element running over all particles </>, 
k, lj, N a , Nb (V is the spatial volume). Standard lep- 
tons and Higgs boson are assumed to be in thermal equi- 
librium while singlet neutrinos have densities given by 
Eq. (|19|l and distribution functions proportional to the 
thermal distribution functions as in Eq. I|l(jfl . The chem- 
ical potentials are equal to zero. The factors 

Jijab = l^{hiahjah* ib h* jb }M a M b (36) 

signal the necessary CP and lepton number violation 
through Yukawa couplings and Majorana masses. Jijab 
are antisymmetric under neutrino flavor exchange a <-> b 
and symmetric under lepton flavor exchange i j. The 
internal product p[ ■ p'j = PiPj — Pi ■ Pj in the integrand 



functions refers to pseudo 4- vectors p[. The p\ space com- 
ponents coincide with the 3-vector linear momentum pi 
and its time component is equal to the momentum ab- 
solute value, pf = pi. The 4- vectors p\ result from the 
spinor wave functions in a thermal bath |55j : ^2 u u =$\. 
In contrast, the time component of the 4-momentumpi is 
the lepton energy, related with the lepton thermal mass 
as usual: p\ = (p 2 + m 2 ) 1 / 2 . 

The generation rate of total standard lepton number 
L = Li is obtained from Eq. H29|) . The constraints 
of Eqs. (|27|1 and (|28|l indicate that only scatterings, or 
better, inverse decays, contribute as L sources: 

(L) s = A 7 (JV a I i -» N b k)s = 

abij 

-5^A 7 (JV„/ j -^0)s. (37) 

aj 

We evaluate the asymmetries including Pauli blocking ef- 
fects and nonzero lepton thermal masses. However, we 
neglect the neutrino masses inside the integrals because 
we are calculating the leading order contributions at rel- 
ativistic temperatures T 3> M a . Neutrino masses appear 
in the constant factors Jijab- All the dependence on the 
neutrino flavors goes in Jij ab and in the neutrino num- 
ber densities n a that are functions of the temperature 
and flavor dependent relaxation temperatures T a . No- 
tice that the following identity holds for the distribution 
functions in Eq. (|3"P 
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UfiQ- + U) = ~rU(l ~ /„° q )(l - fj) ■ (38) 

Under these conditions the generation rate of standard 
total lepton number is 

abij H 

with c ~ 3 (c — 2.6 for a Higgs thermal mass = 
0.6 T and lepton thermal square masses mf — 0.036 T 2 ). 
It is clear that leptogenesis ceases when neutrinos reach 
thermal equilibrium abundances. 



C. Neutrino asymmetries 

In the above equations one sums over both N a neutrino 
helicities. But as emphasized in this paper the leptoge- 
nesis processes generate also asymmetries in the abun- 
dances of the two helicity states and N~ . The helic- 
ity asymmetries are denoted as A a and are defined with 
respect to the cosmological comoving frame. In the fol- 
lowing the rate asymmetries like A 7 (A — > N^Y) denote 
the difference between the rates of the CPT conjugate 
reactions X — > N+Y and X — > N~Y. Wherever the 
neutrino helicity does not appear explicitly a sum over 
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helicities is assumed. The leading order source terms are 
(A b ) s = J2 { A 7^ - N+k) ~ Aj(N+k -> $) 

i 

+ A 7 fo - N+k) - A 7 (N+k - 0)} g 
+ 2 { A 7 (i\y j N+k) - A 7 (iV+? l iVJ,-) (40) 

aij 

+ Aj{N a lj N+k) - A^N+k iV i i )} B . 

When a = 6 in the second summation one gets the correct 
factor of 2 for the helicity flip reactions N^ 1 — ► AT^ . 

The source terms are subject to CPT invariance 
and unitarity conditions namely 

A-y(N+k - 0)s + 53 A 7(A^ -» ^Js = 0, (41a) 

aj 

Aj(N+k ^) s + ^ A 7 (JV+^ -> ATj^s = 0. (41b) 

These constraints eliminate four of the source terms con- 
tained in Eq. (|4T)|) : 

(A 6 ) s = 53 Aj(4> - A^ii) s + A 7 (<£ - iV+r^s 

+ 53 A 7 (A r a ij jv^ii) s + A 7 (iv a r j -> A^+ros- (42) 

From the CPT invariance condition under thermal equi- 
librium and zero chemical potentials, Eq. (|30[) . one de- 
rives also that 

A 7 (^ A^) cq + 53 A 7 (A^ a / j -» JV 6 + Zi)eq = 0, (43a) 

oj 

A 7 (0 AT+y eq + 53 A 7 (JV„r j - -> AT+y eq = 0,(43b) 

aj 

which ensures that the source terms vanish in thermal 
equilibrium. 

The rate asymmetries in Eq. (I29fl contain by definition 
sums over neutrino helicity states and therefore relate to 
the ones of Eqs. (|1*2"|) as 

Aj$^N b k) = Aj(^N+k) 

-A 7 (0 -> N+k) , (44a) 
A 7 (A r a 2j -> JVftZi) = Aj(N a lj -» AT+Zi) 

-Aj(Njj -> N+k) . (44b) 

Replacing this in Eqs. H29|) and (|42|l one obtains the 
leading source terms of the total 'lepton number' i' = 
£L, + £A a : 

(L')s - 253A 7 (0^AT+/ i )s 

bi 

+ 2j2 A 7(N a lj^N+k) s . (45) 

abij 



The asymmetries are calculated from the diagrams of 
Fig. QJ. The A^ positive (negative) helicity states are 
specified with spinors satisfying Eq. ©. The results are 
the following, using the same notations as in Eqs. 1331) - 

A 7 (0 -► N+k) s = 

53 Jij„6 | d$ /^F 3 -F n F t (pa _-^_ m , , (46) 

aj 

A 7 (Ay, - JV+Zi)s = 

Jijafc y d$ fafjF^FiFb (p tt _pt)2_ m a , (47) 

with 

r^^ + ^-^-^Pi,-^. (48) 

Notice that the A^ helicity states as well as the absolute 
3-momenta p,-, Pj, pb in the expressions above are defined 
with respect to the cosmological comoving frame. As be- 
fore we neglect the neutrino masses inside the phase space 
integrals. The integrand function r is finite however, if 
one neglects lepton thermal masses and Pauli blocking 
factors the integration over the fermion angular variables 
yields a null result (recall that the distribution function 
factor in Eq. I|47(l obeys the relation (|38Jl ). This unex- 
pected result means that in that approximation L' is not 
generated at all i.e., the total asymmetry, A = ^^b, 
generated in the neutrino sector cancels exactly the total 
lepton number generated in the standard lepton sector. 
But for nonzero lepton thermal masses the result is finite: 

abij H 

with c ~ 1 (c = 1.2 for Higgs thermal mass = 0.6 T 
and lepton thermal square masses mf = 0.036 T 2 ). 

The comparison with Eq. I|39|) shows that leptogenesis 
generates a total neutrino helicity asymmetry A of oppo- 
site sign and larger than the standard lepton asymmetry 
L so that the source of L' = L + A is of opposite sign 
to L. The neutrino and lepton asymmetries are not sep- 
arately conserved. They are both violated by neutrino 
Yukawa couplings and the lepton number L is violated 
by weak sphalerons as well. These interactions cause the 
interchange of B — V between the singlet neutrino and 
standard lepton and quark sectors but conserve B — L' 
in absence of neutrino Majorana masses. The neutrino 
masses are so responsible for the B — L' dissipation. In 
next section we study this effect and establish the evolu- 
tion of B — L', 



V. WASHOUT PROCESSES 

Weak sphalerons and neutrino Yukawa couplings are 
the only ones that violate the standard partial lepton 
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numbers Li. The processes like <j> — ► iV a ij, iq t — > N a k vi- 
olate the standard total lepton number L — Lj by one 
unit, but this does not mean that they dissipate the total 
lepton number that is generated in leptogenesis. Indeed, 
one has to distinguish between the two neutrino helic- 
ity states N^- The reactions <f> — > N~U, iqt — » N~k, 
4>A^ N~k (A is a gauge boson) -> ^A, -> 
A^+f — > g t Zj, N^q t — > iij violate L by one unit but 
conserve the total 'lepton' number L' = L + A because 
they decrease the total helicity A = ^A a by one unit. 
On the contrary, the same processes with opposite helic- 
ity neutrino states violate L' by two units. They require 
nonzero neutrino Majorana masses and are suppressed 
at relativistic temperatures T 3> M a . The V conserving 
processes dominate the neutrino thermalization and en- 
ter in equilibrium with relaxation temperatures T a given 
by Eq. i|20|) . Below these temperatures the L' conserving 
processes are fast in comparison with the expansion rate 
and thus enforce constraints on the lepton and neutrino 
chemical potentials: 



V4> + Vi ~ Va = 



(50) 



This implies that due to rapid flavor violation all lepton 
doublets U have equal degeneracy parameters rji = rji and 
all positive (negative) helicity neutrinos (N~) have 
equal degeneracy parameters r/ a — rj^ + rji {—rj a ). As a 
result the partial asymmetries Li and A a are not directly 
related with the respective leptogenesis sources but are 
rather determined from the overall B — L' asymmetry by 
the set of chemical potential constraints. 

There are as much constraints as independent fast fla- 
vor changing reactions and this depends on the tempera- 
ture scale [6ll |62| . In the standard model case the latest 
interactions to enter in equilibrium are the Yukawa cou- 
plings and their equilibrium temperatures depend on the 
particular right-handed quark or lepton isosinglet (for a 
detailed discussion see refs. 0, EH)- I* 1 anv case the 
constraints leave Li — L' as a free variable and determine 
the other quantum numbers, in particular Li and A a , as 
proportional to B - L'. In turn, B — L' is completely 
determined by the leptogenesis sources on one hand and 
the dissipation (washout) processes on the other hand. 

B — L' is violated by two units in the AL = 1 , AA = 1 
reactions (f> — > N+li, iqt — > N+li, <f)A — > AT+Zj, and 
crossed channels N~q t — > tZ,, in the AL = 2, A A = 
scatterings <p<f> — > IJj, <f>lj — ► (fili, IJj — > (f><f>, and in 
the AL = 0, AA = 2 scatterings ljj,(j)(t> -> N+N+, 
- <f>N+, l 3 N^ - kN+ 

N-N: 



"b ~ ' lihi4> ( ) ) - All these processes depend on the 
existence of neutrino Majorana masses. In the AL = 2 
case they contribute through the neutrino propagators 
and in the AL = 0, 1 reactions they make possible that 
a neutrino be produced or annihilated with the 'wrong' 
helicity state i.e., helicity opposite to the chirality deter- 
mined by the Yukawa couplings kN a cf). The left-handed 
(right-handed) chiral projection of a positive (negative) 
helicity state with energy E goes as M a /2E in the rela- 
tivistic limit. Thus, the L' violating processes are sup- 



pressed with respect to the respective V conserving chan- 
nels by a ratio going as M^/AE 2 . The thermal aver- 
age of this ratio is about M^/16T 2 for scatterings but 
in the case of cf> — ► N a li decays the center of mass en- 
ergy, rrifp ~ 0.6 T, and average neutrino energy are much 
smaller. The <p ~ * AL"~Z, branching fraction is given by 



IK 



N+k) ^ 8M a 2 



rp2 



(51) 



two orders of magnitude higher than in the case of scat- 
terings. As the AL = 1 scatterings have total rates 
comparable with the Higgs decays, this means that the 
Higgs decays dominate the L' violating rates by two or- 
ders of magnitude over the scatterings. As far as AL = 
and AL = 2 scatterings is concerned, — * N+N^~, 
4>(j) hlj, they are also suppressed by the same 
10~ 2 factor plus an extra factor of \hi a \ 2 with respect to 
the AL' = 2 Higgs decays, because they get one more 
Yukawa coupling than the AL = 1 scatterings. 

From the discussion above one learns that the Higgs 
decays and inverse decays dominate the violation of B — 
L' . Hence, the evolution of B — L' is well described by 
the equation 

B — L' = -(L')s + 2]T{A 7 (7V+Z l - $) 

at 

-A 7 (0 -+ N+k) } , (52) 

where the label '0' indicates that the rates are calcu- 
lated at tree level contrary to the leptogenesis source 
term (L')s given in Eq. I|45|). At temperatures below 
T a neutrinos have thermal equilibrium abundances and 
the transport terms depend on the chemical degeneracy 
parameters as follows: 



B — L' 



-(L')s + ±Y. ^n+iM + m + Va) ■ (53) 



Moreover, the L' conserving neutrino reactions are in 
equilibrium under T a and the constraint of Eq. lfST!|l ap- 
plies, so that rji = rji and r\ a = -q^ + rji. The degen- 
eracy parameters parametrize the particle antiparticle 
number asymmetries. For example, for a ultrarelativis- 
tic fermion / the asymmetry is given in leading order by 
Yf — Yf = 1.83 rjfY eci [6(J,|64(. The constraints imposed 
by the standard model interactions leave one degeneracy 
parameter as free parameter [6ll 162^ . In our case they 
lead to a relation of the form 



B — L' 



-r(Tj<f, + ry ; )Y G( 



(54) 



The precise value of the factor r depends on the set of 
chemical constraints and therefore on the temperature 
scale. It varies from r = 12.8 at temperatures above 
10 12 GeV, where the Yukawa couplings of the isosinglet 
right-handed bu quark and tr lepton are not yet in equi- 
librium, to r = 18.6 under 10 4 GeV, where all quark and 
lepton Yukawa couplings are in equilibrium. 
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The 4> — > N+k partial decay rate is a fraction of the 
4> — » N a li decay rate as indicated in Eq. (|51|1 . On the 
other hand, the Higgs decay rates scale with the first 
power of the temperature and can be related in the same 
way as in Eq. Il2()[| with the Hubble expansion rate H: 

i 

The temperatures T' a are about one half smaller than the 
relaxation temperatures T a because unlike the later T' a 
do not receive contributions from scattering processes. 
Combining everything one rewrites Eq. (|53|) in the form, 

a a 

This shows that when leptogenesis processes are ineffec- 
tive, at temperatures T < T a , B — L' decays exponen- 
tially with T~ 3 and damping constant equal to T^. The 
same happens with the neutrino degeneracy parameters 
Va = V<t> + Vh as determined by Eq. (|54l) . 

It is important to recognize the difference between T* 
and the relaxation temperatures T a . At T a neutrino re- 
actions enter in thermal equilibrium and the lepton and 
neutrino partial quantum numbers Li — B/3 and A a be- 
gin to be rapidly violated. However, it does not mean 
that the standard lepton number L or B — L' start to 
be washed out at T a . It only means that the effective 
couplings liN~(j) enter in equilibrium and the constraints 
rj<t> + Vi ~ Va ~ are enforced. The complete set of chem- 
ical equilibrium constraints force B — L' to distribute 
into B, —L and —A in similar proportions. On the other 
hand the effective couplings liN+4> are n °t yet in equilib- 
rium because the temperatures T a are much larger than 
the neutrino Majorana masses. These couplings enter in 
equilibrium later at the temperature T*. Below that tem- 
perature the constraints r]^ + r/i + i] a = should apply 
on top of the previous constraints + rji — r/ a = 0. It 
means that the degeneracy parameters rj a , r\$ + rji, and 
B — V as well as all asymmetries proportional to B — L' 
are strongly damped below the temperature T+. But not 
above T*. In fact, B — L' is only marginally damped at 
the T a temperature scale. 

For decay constants K a as large as 70, neutrinos enter 
in equilibrium at temperatures T a ss jK a M a ~ 10 M a 
but the B — L' damping constant is two orders of 
magnitude smaller than the T? scale: T? ~ 7777^ T%. 
Moreover, the greater the decay constants K a are the 
smaller is in comparison with T 3 . On the other hand 
T* is quite close to the neutrino mass scale. 

This gives an important lesson. If one ignores neu- 
trino helicity asymmetries, as has ever been done in the 
literature, then one obtains that the temperatures T a set 



the chemical potential constraints rji + 77^ = and define 
the B — L relaxation temperature scale. This is wrong. 
The correct constraints must include the neutrino degen- 
eracy parameters rj a and are given by Eq. I|50ll . Ignoring 
neutrino helicity asymmetries leads to an overestimate 
by two orders of magnitude of the damping rate of the 
lepton and quark asymmetries. 

So far we have considered that singlet neutrinos have 
equilibrium temperatures T a > 10 M a larger than any 
of their masses M a . Another scenario that has been of- 
ten considered for simplicity is the hierarchical scenario 
H II HE EEEl HI where one of the singlet neutrinos, 
Ni, is much lighter than the others. Then, it is argued 
that any asymmetries generated at the heaviest neutrinos 
decaying phases are later washed out by the lightest neu- 
trino AL = 1 reactions before the temperature reaches 
the mass Mi. As we have just shown this is not cor- 
rect because it does not take into account the Ni helicity 
asymmetry. One concludes from Eq. (|57|l that the quan- 
tum number B — L' and all particle asymmetries related 
to it by chemical equilibrium constraints are damped 
with a relaxation temperature given by k, |MfTi 
(in cases where the AL = 2 reactions mediated by off- 
shell neutrinos are not significant). T* is quite close to 
the mass M\ which means that the asymmetries left after 
the heavy neutrino decays remain conserved at tempera- 
tures above the lightest neutrino mass. 

Another possibility is that the lightest neutrino enters 
in equilibrium at a relaxation temperature T\ below the 
other neutrinos masses. Then, the helicity asymmetry 
Ai remains conserved when the heaviest neutrinos decay 
even if the standard lepton asymmetries are completely 
washout by their AL = 1 collisions or mediated AL = 2 
reactions. Later on, the neutrino asymmetry Ai is con- 
verted into a standard lepton asymmetry at tem per atures 
below T\. Preliminary numerical calculations |65j indi- 
cate that as much as 1% of the neutrino asymmetry may 
be converted into a final B — L asymmetry even if there is 
no CP asymmetry associated with the lightest neutrino 
decay. This makes a radical contrast with the traditional 
scenario and will be further investigated. 



VI. CONCLUSIONS 

We studied an aspect of the leptogenesis mechanism 
that has ever been overlooked. We showed that as a 
rule the two helicity states of Majorana neutrinos do not 
have exactly the same abundances, contrary to what has 
been tacitly assumed. The helicity asymmetries defined 
as differences between the two helicity state abundances 
of each neutrino species can be parametrized with appro- 
priate neutrino chemical potentials as any other particles 
asymmetries. The leptogenesis processes generate neu- 
trino helicity asymmetries of the same order of magnitude 
as the standard lepton asymmetries. This is quite natural 
because in the ultrarelativistic limit neutrino Majorana 
masses are negligible and the neutrino - lepton Yukawa 
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couplings conserve a total lepton number assigned as +1 
for right-handed neutrinos and -1 for left-handed neutri- 
nos. 

The neutrino helicity asymmetries participate in the 
system of Boltzmann equations that govern the evolution 
of particle asymmetries and do not decouple because for 
any particular reaction the two neutrino helicity states 
have distinct reaction rates. A reaction where a incom- 
ing left-handed lepton doublet goes with a incoming (out- 
going) neutrino with positive (negative) helicity is sup- 
pressed with respect to the reaction where the same neu- 
trino species is in the opposite negative (positive) helicity 
state. The former is possible only because the neutrino 
has a Majorana mass but its amplitude is suppressed by 
the Lorentz contraction factor. The two reactions rates 
coincide only when the neutrino is at rest. As a result the 
standard lepton and neutrino sectors interchange asym- 
metries with each other which affects significantly the 
transport and dissipation of lepton number. 

It proves convenient to work with the quantum num- 
ber B — L' where L' is the sum of the total standard 
lepton number, L, and the total neutrino helicity asym- 
metry, A (equal to the difference between the total num- 
ber of neutrinos with positive helicity and the total num- 
ber of negative helicity neutrinos). B — L' is conserved 
by sphalerons, its violation by Yukawa couplings is sup- 
pressed by the neutrino mass over temperature ratios, 
and reduces to the standard model B — L number when 
the singlet neutrinos vanish from the Universe. The 
chemical potential constraints enforced by the fast re- 
actions at any given moment set the asymmetries of par- 
ticles in equilibrium as proportional to B — L' . 

We made a detailed analysis of the period when neu- 
trinos are ultrarelativistic which permits simplifying ap- 
proximations. It shows that one has to distinguish be- 
tween the temperature T a at which a neutrino species 
jV a comes into equilibrium and the temperature scale 
where its reactions have a significant damping (washout) 
effect on B — L' and all particle asymmetries propor- 
tional to B — L' (B — L in particular) . A neutrino species 
comes into equilibrium when its AL = 1 reactions be- 
come fast in Hubble rate terms. In the existing literature 
all AL ^ reactions contribute to wash out the B — L 



asymmetry. But this is not right because it does not 
differentiate between the neutrino helicity states. Some 
neutrino helicity configurations conserve L' . These chan- 
nels dominate the rates and exist even in absence of neu- 
trino Majorana masses. The other channels violate L' 
and contribute to wash out B — L' but they are subdom- 
inant and are suppressed by the second power of neutrino 
Lorentz contraction factors. A consequence of this is that 
the contribution of AL = 1 reactions to the B — L 1 damp- 
ing rate is two orders of magnitude smaller than expected 
if the neutrino helicity asymmetries are ignored. 

In section IIVI we studied the generation processes in 
the period of neutrino thermal production. It turns out 
that the B — L 1 asymmetry generated directly in the neu- 
trino helicity sector is of opposite sign and larger than 
the asymmetry generated in the standard lepton sector. 
We did not attempt to study the decaying phase of sin- 
glet neutrino(s). The calculations are more involved be- 
cause the temperature is of the order of the neutrino 
mass(es). But one can identify one important difference 
with the traditional hierarchical scenario, where one of 
the singlet neutrinos is much lighter than the others: the 
asymmetries left when the heaviest neutrinos vanish from 
the Universe are not necessarily washed out by the fast 
AL = 1 reactions of the lightest neutrino before the tem- 
perature reaches its mass scale. Moreover, even if the 
heaviest neutrino AL = 1 collisions or AL = 2 medi- 
ated reactions washout the standard lepton and Higgs 
boson asymmetries, the helicity asymmetry carried by 
the lightest neutrino may survive if it enters in equilib- 
rium at a relaxation temperature below the heaviest neu- 
trino masses. The neutrino asymmetry is later converted 
into a standard global B — L asymmetry at temperatures 
close to its mass. This has been supported by numerical 
calculations and will be reported elsewhere |65j . 
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